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(54) Title: METHOD, DEVICE AND SOFTWARE PACKAGE FOR SMOOTHING A SUBSURFACE PROPERTY 

(54) Titre : PROCEDE, DISPOSITIF ET PRODUTT -PROGRAMME DE LISSAGE D'UNE PROPRIETE DE SUBSURFACE 

(57) Abstract: The invention concerns a method for 
smoothing a subsurface property in a geological structure 
represented by seismic measurements comprising a step 
(100) which consists in digital modelling by continuous local 
seismic traces, calculating an optimal offset and defining a 
conditional neighbourhood of a reference central continuous 
local seismic trace; a step (101) which consists in selecting 
the property to be smoothed on a conditional neighbourhood; 
a third step (102) which consists in substituting properties 
of the conditional neighbourhood with offset properties; and 
a fourth step (103) which consists in selecting an average of 
the properties offset on the conditional neighbourhood at step 
(102). 
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(57) Abrege : Un procSde' de lissage d'une propri&e* de sub- 
surface dans une structure geblogique representee par mesures 
sismiques comprend une 6 tape 100 de modelisation numerique 
par traces sismiques locales continues, de calcul du dScalage 
optimal et de definition du voisinage condiuonnel d'une trace 
sismique locale continue "centrale" de reference ; une 6tape 
101 de choix de la proprigte* a lisser sur un voisinage condiuon- 
nel ; une troisieme 6tape 102 de remplacement des propri6t6s 
du voisinage conditionnel par des proprtetes decalees ; et une 
quatrieme 6tape 103 de choix d'une moyenne des propri6t6s 
d6calees sur le voisinage conditionnel a T6tape 102. 
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Vmr*A6. dispositif et produ ct ™ pramine de lissage d'une propriete de 
subsurface . 

L'invention est relative a un proced6 de Ussage d'une propriete de subsurface 
dans une structure geologique representee par des mesures sismiques. 
5 L'invention est egalement relative a un dispositif de lissage d'une propriete de 

subsurface dans une structure geologique representee par des mesures sismiques. 

L'invention est enfin relative a un produit-programme d'ordinateur permettant 
le fonctionnement d'un dispositif programmable de lissage d'une propriete de 
subsurface dans une structure geologique representee par des mesures sismiques. 
10 Le document WO 01/63323 Al decrit un procede de traitement de donnees 

sismiques comprenant les etapes consistant a : obtenir un volume de donnees 
sismiques couvrant un volume pred&ermin6 de terre ; determiner pour chaque voxel 
du volume de donnees sismiques l'orientation locale des donnees sismiques ; 
d6terminer pour chaque voxel s'il existe un bord dans le voisinage, et executer une 
15 operation de Ussage sur chaque voxel dans le volume de donnees sismiques, dans 
lequel la direction de l'operation de Ussage est l'orientation locale des donnees, et dans 
lequel l'operation de Ussage ne depasse pas le bord, de maniere a obtenir un volume de 
donnees sismiques traitees, ou la quantit6 associee a chaque voxel dans le volume de 
donnes traitees est le resultat obtenu par l'execution de l'operation de Ussage dans le 
20 voxel dans le volume de donnees sismiques. 

Le document WO 02/13139 Al decrit un procede de traitement d'images 
sismiques comprenant les etapes : obtenir un ensemble de donnees d'une image initiale 
bidimensionnelle ou tridimensionneUe, dans lequel chaque element de l'ensemble de 
donnees est l'intensite initiale d'image du point de l'image ; calculer pour chaque point 
25 les derivees partielles de l'element dans n directions pour obtenir un ensemble de 
valeurs derive des derivees partieUes; calculer pour chaque point une matrice carree 
structural symetrique a partir des valeurs des derivees partieUes; effectuer une 
iteration point a point ponder6e par une variable voisine de zero lorsque le point est 
voisin d'un bord et voisine de 1 lorsqu'U est tres eloigne d'un bord ; et repeter ces 
30 etapes un nombre de fois pour obtenir l'image traitee. 

Un premier but de l'invention est d'ameliorer le Ussage des proprietes de 
subsurface dans une structure geologique representee par des mesures sismiques. 
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Un deuxi&ne but de 1'invention est de permettre un lissage simple et rapide des 
proprietes de subsurface, tout en ne lissant pas les discontinuites. 

L'invention a pour objet un precede de lissage d'une propri6te de subsurface 
dans une structure geologique representee par mesures sismiques, dans lequel on 
5 construit une fonction continue Si j)k (t) par interpolation ou approximation des traces 
sismiques discretes d'une matrice sismique multidimensionnelle, ladite fonction etant 
designee comme "trace sismique locale continue", comportant les etapes suivantes : 

a) . uuliser comme decalage optimal h. ii>pq ^ de deux traces sismiques locales 
continues voisines S ij)k (t) et S^t), la valeur de d6calage rendant maximale leur 

1 0 fonction de correlation ; 

b) retenir comme voisinage conditionnel d'une trace sismique locale continue 
"centrale" S ii)k (t) le sous-voisinage consistant en des traces voisines Spq^t) 
correspondant a des d6calages optimaux UimM associes a des correlations R&p.^ij.pq*) 
superieures a un seuil predetermine compris entre 0 et 1 ; 

15 c ). choisir une propri6te de subsurface a lisser au voisinage conditionnel d'un 

point (i, j, k) d'une trace sismique locale continue "centrale" de reference ; 

d) . decaler les proprietes de subsurface du voisinage conditionnel en translatant 
la variable courante de la valeur du d6calage optimal hij >pq ,k; 

e) . prendre comme valeur liss6e au point (i, j, k) une moyenne des proprietes de 
20 subsurface d6cal6es a l'&ape d). 

Selon d'autres caracteristiques de l'invention : 

- la moyenne des proprietes de surface de l'etape e) est une moyenne qui peut 
Stre ponderee, par exemple par la valeur de correlation maximale correspondant au 
d£calage optimal, 

25 - la moyenne des proprietes de subsurface de l'etape e) est choisie dans 

l'ensemble suivant : moyenne arithmetique, moyenne geometrique, moyenne 
harmonique, 

- on choisit comme propriete de subsurface a lisser ramplitude reflechie et 
captee par des geophones, 
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- si nScessaire, V ensemble du processus de lissage est r€p&6 un certain nombre 
de fois pour ameliorer le lissage, 

. on visualise des coupes de matrice multidimensionnelle de propri6te lissee sur 
un ecran de visualisation. 

5 L ! invention a egalement pour objet un dispositif pour la mise en ceuvre du 

procedS selon l'invention, comportant des moyens pour utiliser comme decalage 
optimal de deux traces sismiques locales continues voisines la valeur de d6calage 
rendant maximale leur fonction de correlation, des moyens pour dScaler les propri6t6s 
de subsurface du voisinage conditionnel en translatant la variable courante de la valeur 

10 du decalage optimal hy,^ des moyens pour choisir une propriete de subsurface k 
Usser au voisinage conditionnel d'un point (i, j, k) d'une trace sismique locale continue 
"centrale" de reference, des moyens pour dealer les propri&es de subsurface du 
voisinage conditionnel en translatant la variable courante de la valeur du decalage 
optimal hij,pq,k et des moyens pour prendre comme valeur lissee au point (i, j, k) une 

1 5 moyenne des propri&es de subsurface d6cal6es k l'etape d). 

Selon d'autres caractfristiques de l'invention : 

-le dispositif comporte des moyens de m6morisation et des moyens de 
visualisation de parametres sismiques d6termin6s a l'aide du proc&te selon l'invention. 

L'invention a enfin pour objet un produit-programme d'ordinateur, comportant 
20 des Elements de code de programme pour executer les 6tapes d'un proc6de selon 
Tinvention, lorsque ledit programme est execute par un ordinateur. 

L'invention sera mieux comprise grace k la description qui va suivre donnee k 
titre d'exemple non limitatif en reference aux dessins annexes dans lesquels : 

_ La figure 1 repr6sente sch&natiquement une matrice sismique 
25 tridimensionnelle et une trace sismique locale continue. 

- La figure 2 repr6sente schematiquement, un exemple de voisinage local N ijf k 
consistant en un ensemble de trace sismique locale continues S pqfk (t) voisines de la 
trace centrale de r6ference S ijj£ (t) elle m&ne situ6e a la verticale du geophone G {j et 
interpolant les amplitudes sismiques autour de t=k. 
30 - La figure 3 repr6sente sch&nafiquement, une section verticale de voisinage 

local conditionnel Nij, k (r) d'une trace sismique locale continue de r6f6rence S ij>k (t). 
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- La figure 4 represente schematiquement, une section verticale d'une matrice 
sismique tridimensionnelle avec un decalage optimal h ijjPq ^ et un voisinage 
conditionnel Nij^(r). 

- La figure 5 represente schematiquement, une section verticale d'une matrice 
5 sismique tridimensionnelle analogue a la figure 4 avec presence d'une faiUe ou 

discontinue 

La figure 6 represente schematiquement, un organigramme fonctionnel d'un 
procede selon 1'invention. 

En reference aux figures 1 a 5, les el&nents identiques ou fonctionnellement 
1 0 Equivalents sont d6sign6s ou reperes de maniere identique. 

Sur la figure 1, une matrice sismique tridimensionnelle est obtenue par releve 
de mesures enregistrees par des geophones G i} disposes sur un reseau x,y en des points 
de coordonnees i, j. Le relev6 des mesures echantillonnees dans le temps est repr6sente 
suivant un axe t descendant repr6sentatif de la profondeur ou d'une verticale 
15 descendante & partir de la surface de la terre ou de la mer. Les mesures sont 
caract6ris6es par leur amplitude, par exemple une amplitude relevee par le gSophone 
Gg au temps ou k la profondeur d'echantillonnage t k . La mesure discrete effectu6e par 
le g6ophone Gi, au temps ou & la profondeur t k est appelee amplitude sismique S ij>k . 

L'ensemble des amplitudes sismiques correspondant k un g6ophone Gy de 
20 coordonnees i, j est une matrice unidimensionnelle (Sgi, Sij 29 — » S ijk ,. . .., Sijn) appeiees 
trace sismique discrete, car cette matrice unidimensionnelle correspond k la trace selon 
le point de coordonnees horizontals i, j de la matrice sismique tridimensionnelle 
obtenue par mesures sismiques. 

L'axe vertical t orients selon une verticale descendante d6signe usuellement le 
25 temps, mais peut 6galement 6tre traite pour repr^senter une profondeur k partir de la 
surface. 

L'invention conceme aussi bien l'application k une troisteme coordonn6e t 
representative du temps, qu'&. une troisieme coordonn6e t representative de la 
profondeur. 

30 A partir de la trace sismique discrete situ6e k la verticale d'un geophone Gg on 

d6finit, par interpolation ou approximation des valeurs discretes autour de t=tk=k, une 
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fonction continue S ij)k (t) qui est designee comme "trace sismique locale continue". Les 
methodes d'approximation ou ^interpolation de valeurs discretes pour engendrer une 
fonction continue sont nombreuses, et comprennent notamment les interpolations ou 
approximations polynomiales, ainsi que les interpolations ou approximations 
5 trigonometriques polynomiales. 

Toute autre variante ^interpolation ou d'approximation fournissant une 
fonction continue peut egalement etre appliquee a la presente invention pour fournir 
une "trace sismique locale continue". 

Sur la figure 2, plusieurs traces sismiques locales continues definissent un 
10 exemple de voisinage d'une trace sismique locale continue "centrale" Syj^t) de 
reference. Le voisinage d'une trace sismique locale continue S ij>k (t) est ddfini comme 
l'ensemble des traces sismiques locales continues dont les indices spatiaux horizontaux 
p,q sont voisins des indices spatiaux horizontaux i, j de la trace sismique locale 
continue de reference. 

15 A titre d'exemple, les coordonnees spatiales horizontales p,q correspondant au 

geophone G pq sont voisines des coordonnees spatiales horizontales i, j correspondant 
au geophone Gy si les valeurs absolues des differences i-p et j-q sont inferieures a des 
entiers donnes, par exemple a 2. 

Dans ce cas, comme represents a la figure 2, la trace sismique locale continue 
20 S ij>k (t) est associee a huit traces sismiques locales continues voisines entourant la trace 
sismique locale continue "centrale" Sij,k(t). 

Dans le cas de traces sismiques locales continues produites par des mesures 
sismiques, le profil des horizons geologiques introduit des decalages verticaux entre 
des traces sismiques locales continues voisines. Dans le but de determiner des relations 
25 entre deux traces sismiques continues locales voisines S ijJc (t) et Spq,k(t) centrees sur la 
m&ne coordonn6e verticale d'6chantillonnage t=t k =k et correspondant a des 
coordonnees spatiales i, j et p,q differentes, on calcule la fonction de correlation 
Rij,pq*(h) des deux traces sismiques locales continues. 

La fonction de correlation de deux traces sismiques locales continues voisines 
30 est obtenue par la formule suivante 
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oix le num6rateur correspond a la fonction de covariance de Sij^t) et Spq^t) 

tk+A 

obtenue par l'expression suivante Qj,pq*(h) = J S ijk (t) . (t + h) . dt . 



Dans cette integrate dSfinissant Qj^Oi) , le param&re A d6finit vine « fenetre verticale 
d'investigation » autour de t=tk=k 

Par exemple, si les traces sismiques locales continues Sy,k(t) et Spq^(t) sont des 
polyndmes trigonometriques de la forme suivante interpolant les donnes sismiques 



m 

Sij,k(0= £af .cos(sa>t) + bf .sin(s^t) 

m 

io Sij^t) = £af .cos(s<»t)+bf .sin(sa>t) 



alors on peut montrer que la fonction de covariance Cij (Pq>k (h) est elle meme un 
polynSme trigonom&rique de la fonne suivante ou les coefficients A s lipq,k et B s 8 ' pq ' k 
dependent des coefficients a,**, a*™*, b a U et b™* : 

15 C ij>pq)k (h) = jAf^-cosCso^ + Br-'-sinCs^h) 

Un r6sultat math&natique classique et connu est que la fonction de correlation 

R.. ,(h\ = C ij>pq^ ^ traduit une similarit6 des fonctions corr616es 

lorsque cette fonction de correlation se rapproche de 1. 

20 

L'etude des fonctions de correlation des traces sismiques locales continues 
Si j>k (t) et S pq>k (t) permet de definir comme d6calage optimal la valeur hi itPq>k de h 
correspondant au maximum de la fonction de correlation Rij,pq,k(b) , c'est-a-dire au 
maximvnn de la fonction de correlation le plus proche de 1. Contrairement a l'art 
25 existant base sur une formulation discrete de la fonction de correlation Ri,- >pqk k(h), 
l'utilisation d'une formulation continue des traces sismiques et done de la fonction de 
correlation permet d'obtenir un decalage optimum qui n'est pas contraint a etre un 
multiple entier du pas d'echantillonnage suivant l'axe vertical correspondant a la 
variable t. 
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La definition du decalage optimal hij >pq> k de deux traces sismiques locales 
continues voisines pennet d'obtenir une premiere approximation des horizons passant 
par cette trace sismique locale continue de reference Sij^(t), comme repr6sente aux 
figures 3 a 5. Contrairement k Tart existant, le fait d'utiliser des decalages optimaux 
5 qui ne sont pas contraint a etre un multiple entier du pas d^chantillonnage suivant 
l'axe vertical pennet d'6viter des erreurs connues sous le nom « d 5 aliasing ». 

Sur ces figures 3 a 5, on definit un voisinage conditionnel N^r) de la trace 
sismique continue locale de reference Sij^(t) qui est un sous-voisinage du voisinage 
initial des traces sismiques locales continues retenues pour le calcul de correlation et 
10 de definition des decalages optimaux. Le voisinage conditionnel Nij^(r) est choisi de 
mani£re que, pour toute trace sismique continue locale Spq^t) appartenant a N^r), le 
decalage optimal hij, pq jc correspondant est tel que la correlation Rij,pq,k(hij,pq,k) entre 
Sij,k(t) et S pquk (t) est sup6rieure a un seuil predetermine r compris entre 0 et 1. 

Comme on le voit sur la figure 5, la fixation du seuil r pennet egalement 
15 d'englober dans un voisinage conditionnel une discontinuite ou une faille geologique, 
ce qui constitue un avantage important par rapport a Tart anterieur. 

Les dispositions precit6es fournissent ainsi une modeiisation continue 
permettant la mise en oeuvre de Pinvention. 

En reference k la figure 6, un proc6de de lissage ou propriete de subsurface 
20 dans une structure geologique repr6sente par des mesures sismiques comporte une 
premiere etape 100 de modeiisation numerique, comme decrit en reference aux figures 
14 5, pour definir des traces sismiques locales continues d'une matrice sismique 
multidimensionnelle, calculer les decalages optimaux de traces sismiques locales 
continues voisines et definir les voisinages conditionnels des traces sismiques locales 
25 continues "centrales" de reference. 

Apres cette premiere etape 100 de modeiisation numerique, une 6tape 101 
pennet de choisir la propriete k lisser sur un voisinage conditionnel. Cette propriete k 
lisser sur un voisinage conditionnel peut consister en toute fonction 
multidimensionnelle de la subsurface h analyser. 

30 En particulier, la propriete k lisser peut simplement consister en Pamplitude du 

signal sismique reflechi par les horizons de la subsurface en direction des geophones 
de mesure. 
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A l'6tape 102, on remplace les propriet6s du voisinage conditionnel par les 
propri6t6s d6cal6es de la valeur du decalage optimal hij jPq ^ de maniere a constituer un 
ensemble de proprietes relatives au meme horizon. 

Ainsi, dans le cas d'une trace locale sismique continue centrale de reference 
5 Sij^(t), on remplace toutes les valeurs Spq^(t) du voisinage conditionnel par les valeurs 

Spq^(t+hij iPq ,k) • 

On effectue ensuite k T6tape 103 une moyenne des valeurs decalees 
S pq ,k(t+hij,pq*) etablies a T6tape 102 et on choisit comme valeur de la propriete au point 
de coordonn6es spatiales i j et de coordonnee temporelle ou de profondeur t=tk=k , la 
1 0 moyenne ainsi calculee. 

L'invention s'applique k plusieurs variantes d'etablissement de moyennes : on 
peut par exemple utiliser une moyenne ponderee ou non ponder6e; ou, 
alternativement, on peut egalement utiliser une moyenne arithmetique, une moyenne 
geometrique, une moyenne harmonique , ponderee ou non ponder6e. 

1 5 Dans le cas d'une moyenne pond6r6e, ont peut avantageusement utiliser comme 

facteur de pond6ration la valeur maximale de la fonction de correlation 

Rij >pq * (hij,pq,k) y en refl6tant ainsi le degr6 de confiance accordable k chaque 
decalage optimal correspondant k xin maximum d'une fonction de correlation. 



20 permet ainsi d'ameiiorer le contraste et la precision des propri6t6s de subsurface dans 
une matrice sismique multidimensionnelle. 

Si n6cessaire 5 afin d'ameiiorer le lissage, l'ensemble du processus decrit par la 
presente invention peut etre repute sur le lissage obtenu lors d'une iteration precedente. 
L'etape 104 constitue une etape iterative dans laquelle l'operateur peut dSfinir le 
25 nombre d'iterations adequates k la resolution de son probl&me. 

L'invention est de preference mise en oeuvre sur un dispositif comportant des 
moyens pour accomplir successivement les etapes du proc6d6 d6crit en reference k la 
figure 6. 



30 parametres sismiques ou propri6t6s obtenues k retape 105 et des moyens de 
memorisation de programme et de calculs intermediates. 



Le proc6de de lissage ou de filtrage d'une propriete selon la presente invention 



En particulier, ce dispositif comporte des moyens de visualisation des 
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De preference, ce dispositif est un dispositif programmable commande par un 
produit-programme d'ordinateur, realise pour mettre en oeuvre un procede selon 
Tinvention. 

L'invention d6crite en reference a plusieurs objets particuliers n'y est millement 
limit6e mais couvre au contraire toute modification de forme et toute variante de 
realisation dans le cadre et Tesprit de Tinvention, Tessentiel 6tant dhitiliser des valeurs 
de propriety recalees par rapport aux decalages optimaux determines comme valeurs 
correspondant a la maximisation des fonctions de correlation de traces sismiques 
locales continues ou meme discretes. 
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R RVENDIC ATIONS 

1. Proc6de de lissage d'une propriete de subsurface dans une structure 
g6ologique representee par mesures sismiques, dans lequel on construit une fonction 
continue S ij>k (t) par interpolation ou approximation des traces sismiques discretes d'une 
5 matrice sismique multidimensionnelle, ladite fonction Sy^Ct) etant d6sign6e comme 
"trace sismique locale continue", comportant les etapes suivantes : 

a) , utiliser comme d6calage optimal de deux traces sismiques locales continues 
voisines S^t) et Sp^t), la valeur de d6calage rendant maximale leur fonction de 
correlation, ce decalage optimal n' etant pas contraint a etre un multiple entier du pas 

10 d' 6chantillonnage vertical ; 

b) . retenir comme voisinage conditionnel d'une trace sismique locale continue 
"centrale" Sij^(t) le sous-voisinage consistant en des traces voisines S pq ^(t) 
correspondant k des decalages optimaux associ6s k des correlations Rij sP q*(h) 
superieures k un seuil predetermine compris entre 0 et 1 ; 

15 c ). choisir une propri6t6 de subsurface k lisser au voisinage conditionnel d'un 

point (i, j, k) d'une trace sismique locale continue "centrale" de reference ; 

d) . d6caler les propri6t6s de subsurface du voisinage conditionnel en translatant 
la variable courante de la valeur du decalage optimal (hij^q^ ; 

e) . prendre comme valeur liss6e au point (i, j, k) une moyenne des propriet6s de 
20 subsurface d6cal6es a l'etape d). 

2. Proced6 selon la revendication 1, dans lequel la moyenne des propri6t6s de 
surface de l'etape e) est une moyenne pond6r6e, par exemple par la valeur de 
correlation maximale correspondant au decalage optimal. 

3. Proced6 selon la revendication 1 ou 2, dans lequel la moyenne des propri£t6s 
25 de subsurface de l'etape e) est choisie dans 1'ensemble suivant : moyenne arithm6tique, 

moyenne g6om£trique, moyenne harmonique ponder6e ou non. 

4. Proced6 selon la revendication 1, dans lequel on choisit comme propriety de 
subsurface k lisser, l'amplitude r6fl6chie et capt6e par des geophones. 

5. Proc6de selon rune quelconque des revendications precedentes, dans lequel 
30 on itdre le processus de lissage en Pappliquant sur le resultat d'un lissage anterieur. 



WO 2004/001451 





7FR2003/001755 



11 



6. Proc6de selon la revendication 1 ou la revendication 5, dans lequel on 
visualise une matrice multidimensionnelle de propriety liss6e sur tin ecran de 
visualisation. 

7. Dispositif pour la mise en oeuvre du proc6de selon Tune quelconque des 
5 revendications 1 a 4, comportant des moyens pour utiliser comme d6calage optimal de 

deux traces sismiques locales continues voisines la valeur hj t pq^ de decalage rendant 
maximale leur fonction de correlation, des moyens pour decaler les propriety de 
subsurface du voisinage conditionnel en translatant la variable courante de la valeur du 
decalage optimal hij jPq> k , des moyens pour choisir une propriety de subsurface k lisser 
10 au voisinage conditionnel d'un point (i, j, k) d'une trace sismique locale continue 
"centrale" de reference, des moyens pour d6caler les propriety de subsurface du 
voisinage conditionnel en translatant la variable courante de la valeur du decalage 
optimal hij >pqj k et des moyens pour prendre comme valeur liss6e au point (i, j 5 k) une 
moyenne des propri6t6s de subsurface d6cal6es h l'6tape d). 

15 8. Dispositif selon la revendication 7, comportant des moyens de memorisation 

et des moyens de visualisation de parametres sismiques d6termin6s h Taide du proced6 
selon Tune quelconque des revendications 1 & 6. 

9. Produit-programme d'ordinateur, comportant des Elements de code de 
programme pour ex6cuter les etapes du procede selon Tune quelconque des 
20 revendications 1 & 5, lorsque ledit programme est ex6cut6 par un ordinateur. 



10. Produit-programme d'ordinateur, comportant des elements de code de 
programme pour ex6cuter les 6tapes du proc6d6 selon la revendication 6, lorsque ledit 
programme est execute par un ordinateur. 



INTERNATION 



RCH REPORT 



IntematlJ^^^pli cation No 

PCT/fHB/01755 



A. CLASSIFICATION OF SUBJECT MATTER . _ 

IPC 7 GOlVl/30 G01V1/28 



According to International Patent Ctassiftcallon (IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

IPC 7 G01V 



Documentation searched other than minimum documentation to the extent that such documents are Included in the fields searched 



Electronic data base consulted during the international search (name of data base and, where practical, search terms used) 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category c 



Citation of document, with Indication, where appropriate, of the relevant passages 



Relevant to daim No. 



US 5 615 171 A (HILDEBRAND HAROLD A) 
25 March 1997 (1997-03-25) 
column 5, line 31 -column 11, line 45 
figures 1-10 

US 6 138 076 A (GRAF KERMIT E ET AL) 
24 October 2000 (2000-10-24) 
column 5, line 4 -column 6, line 67 
figures 19,20,23-28 

US 4 672 546 A (FLINCHBAUGH BRUCE E) 
9 June 1987 (1987-06-09) 
column 2, line 33 -column 8, line 17 
figures 1-10 

-/-- 



1-10 



1-10 



1-10 



Further documents are listed in the continuation of box C. 



ID 



Patent family members are listed in annex. 



0 Special categories of cited documents : 

■A" document defining the general state of the art which Is not 

considered to be of particular relevance 
"E" earlier document but published on or after the international 

filing date 

•L" document which may throw doubts on priority ciaim(s) or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 

"O" document referring to an oral disclosure, use, exhibition or 
other means 

"P" document published prior to the international filing date but 
later than the priority date claimed 



r document pubOshed after the international fifing date 
or priorfty date and not In conflict with the application but 
cited to understand the principle or theory underlying the 
Invention 

■X" document of particular relevance; the claimed Invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document Is taken alone 

■Y" document of particular relevance; the claimed invention 

cannot be considered to involve an inventive step when the 
document Is combined with one or more other such docu- 
ments, such combination being obvious to a person sWDed 
In the art 

•&' document member of the same patent family 



Date of the actual completion of the International search 



6 November 2003 



Date of mailing of the International search report 



13/11/2003 



Name and mailing address of the ISA 

European Patent Office, P.B. 5818 Patentlaan 2 
NL - 2280 HV Rijswqk 
TeL (+31-70) 340-2040. Tx. 31 651 epo nl, 
Fax: (+31-70) 340-3016 



Authorized officer 



Modesto, C 



Form PCT/ISA/210 (second sheet) (July 1992) 



INTERNATIOI 



EARCH REPORT 



lntematl<^^^pIlcatton No 

PCT/fPf3/01755 



C(Contlnuatlon) DOCUMENTS CONSIDERED TO BE RELEVANT- . 



Category * Citation of document, with indication .where appropriate, of the relevant passages 



Relevant to claim No. 



US 5 930 730 A (KIRLIN R LYNN ET AL) 
27 July 1999 (1999-07-27) 
column 5 -column 8 
figures 1-10 

US 2002/022930 Al (DALLEY RICHARD MAPES 
ET AL) 21 February 2002 (2002-02-21) 
cited in the application 
page 1, paragraph 16 -page 4, paragraph 59 
figures 1-6 



1-10 



1-10 



US 6 151 555 A (PEPPER RANDOLPH E F 
AL) 21 November 2000 (2000-11-21) 
column 3, line 45 -column 5, line 38 
figures 1-39 



ET 



1-10 



Form PCT7ISA/210 (continuation of second sheet) (July 1692) 



INTERN ATI 

Inf. 



JATI^A 

orma^^Hi 



L SEARCH REPORT 

patent family members 



Intematit 

PCT/I 



illcatlon No 

1/01755 



Patent document 




Publication 




• Pfltont family 




Publication 


cited In search report 




date 




member(s) 




date 


US 5615171 


A 


25-03-1997 


US 


5432751 


A 


11-07-1995 






US 


5251184 


A 


05-10-1993 






US 


5153858 


A 


06-10-1992 






AT 


129577 


T 


15-11-1995 






AT 


190136 


T 


15-03-2000 






CA 


2088501 


Al 


10-01-1993 






CN 


1070267 


A ,B 


24-03-1993 






DE 


69205673 


Dl 


30-11-1995 






DE 


69205673 


T2 


20-06-1996 






DE 


69230731 


Dl 


06-04-2000 






DE 


69230731 


T2 


03-08-2000 






DK 


548341 


T3 


04-12-1995 






DK 


652447 


T3 


24-07-2000 






EP 


0548341 


Al 


30-06-1993 






FP 

Lr 


0652447 


A2 


10-05-1995 






ES 


2144480 


T3 


16-06-2000 






GR 


3033559 


T3 


29-09-2000 






NO 


930843 


A 


07-05-1993 






RU 


2107931 


CI 


27-03-1998 






UO 


9301508 


Al 


21-01-1993 


Uo DIjoU/D 


A 


24-10-2000 


LIS 


6014343 A 


11-01-2000 






CA 


2219216 Al 


30-04-1998 






FR 


2755244 Al 


30-04-1998 






FR 


2760275 Al 


04-09-1998 






fiR 

VID 


2320968 A ,B 


08-07-1998 






NO 


975004 A 


04-05-1998 


US 4672546 


A 


09-06-1987 


CA 


1249876 Al 


07-02-1989 






EP 


0181216 A2 


14-05-1986 






NO 


854447 A 


09-05-1986 


US 5930730 


A 


27-07-1999 


US 


5563949 


A 


08-10-1996 






AU 


709621 


B2 


02-09-1999 






AU 


7382996 


A 


28-04-1997 






CA 


2204168 


Al 


10-04-1997 






CN 


1166207 


A 


26-11-1997 






EP 


0796442 


Al 


24-09-1997 






NO 


971801 


A 


05-06-1997 






UO 


9713166 


Al 


10-04-1997 






US 


RE38229 


El 


19-08-2003 






AU 


696742 


B2 


17-09-1998 






AU 


4133396 


A 


03-07-1996 






CA 


2179901 


Al 


20-06-1996 






CN 


1138902 


A ,B 


25-12-1996 






EG 


20609 


A 


30-09-1999 






EP 


0736185 


Al 


09-10-1996 






NO 


962731 


A 


11-10-1996 






RU 


2144683 


CI 


20-01-2000 






UO 


9618915 


Al 


20-06-1996 






US 


5838564 


A 


17-11-1998 






RU 


2187130 


C2 


10-08-2002 


US 2002022930 


Al 21-02-2002 


AU 


5616501 


A 


03-09-2001 






CA 


2398854 Al 


30-08-2001 






UO 


0163323 Al 


30-08-2001 






EP 


1257850 Al 


20-11-2002 






OP 


2003524192 T 


12-08-2003 



Form PCT/ISA/210 (patent family annex) (July 1992) 



INTERNAT 



RNAT^WAL S 
Informa^^Bpatent 



AL SEARCH REPORT 

family members 



Heation No 



Internatlc 

PCT/FrPR/01755 



Patent document 
cited in search report 



Publication 
date 



Patent family 
member(s) 



Publication 
date 



US 2002022930 


Al 

#1 J. 


NO 


20024030 A 


23-08-2002 


US 6151555 


A 


21-11-2000 AU 


3866300 A 


28-09-2000 






CA 


2362848 Al 


14-09-2000 






EP 


1198776 Al 


24-04-2002 






GB 


2370878 A ,B 


10-07-2002 






NO 


20014353 A 


07-09-2001 






WO 


0054207 Al 


14-09-2000 



Form PCT/1SA/210 (patent family annex) (July 1998) 



RAPPORT DE RE' 



CHE INTERNATIONALE 



Demand 



■ IF 



itlonale No 



PCT/RW3/01755 



A. CLASSEMENT DE L'OBJET DE LA DEMANDE 

CIB 7 G01V1/30 G01V1/28 



Seton la classification Internationale des brevets (CIB) ou a la tote selon la classification naUonate et la CIB 



B. DOMAINES SUR LESQUELS LA RECHERCHE A PORTE 



Documentation mlnimale consultee (systeme de classification sulvi des symboles de dassement) 

CIB 7 G01V 



Documentation consulted autre que la documentation mlnimale dans la mesure ou ces documents retevent des domaines sur lesquels a ports la recherche 



Base de donnees electronlque consultee au cours de la recherche Internationale (nom de la base de donnees, et si realisable, termes de recherche utilises) 



C. DOCUMENTS CONSIDERES COMME PERTINENTS 



Categorie • Identification des documents cites, avec, la ces echeant, rindication des passages pertinents 



no. des revendicaHons vlsees 



US 5 615 171 A (HILDEBRAND HAROLD A) 
25 mars 1997 (1997-03-25) 
colonne 5, Hgne 31 -colonne 11, Hgne 45 
figures 1-10 

US 6 138 076 A (GRAF KERMIT E ET AL) 
24 octobre 2000 (2000-10-24) 
colonne 5, Hgne 4 -colonne 6, Hgne 67 
figures 19,20,23-28 

US 4 672 546 A (FLINCHBAUGH BRUCE E) 
9 ju1n 1987 (1987-06-09) 
colonne 2, ligne 33 -colonne 8, Hgne 17 
figures 1-10 

-/- 



1-10 



1-10 



1-10 



[3 v <* la 



suite du cadre C pour la fin de la Qste des documents 



Les documents de families de brevets sont Indlques en annexe 



° Categories speciales de documents cites: 

"A" document definissant I'etat general de la technique, non 

considers comme partfcutlerement pertinent 
'E" document anterieur, mais publie a la date de depot international 

ou apres cette date 
"L" document pouvant jeter un doute sur une revendication de 

priortte ou cite pour determiner la date de publication cTune 

autre citation ou pour une raison spedale (telle qu'indique*e) 
'O* document se r6ferant a une divulgation orale, a un usage, a 

une exposition ou tous autres moyens 
"P* document publie avant la date de d6p6t international, mais 

posterieurement a la date de priority revendlquee 



T document ulterleur publie apres la date de dep6t international ou la 
date de priortte et n'appartenenant pas a fetal de la 
technique pertinent, mais cite pour comprendre le principe 
ou la theorfe constfluant la base de Invention 

■X s document partlculierement pertinent; rinven tion revendlquee ne peut 
etre consideree comme nouvetle ou comme impiiquant une activfle 
inventive par rapport au document considers isolement 

■Y' document particullerement pertinent; rinven tion revendlquee 
ne peut 6tre consideree comme impiiquant une activlte Inventive 
lorsque le document est associe a un ou plusieurs autres 
documents de meme nature, cette combinalson etant evidente 
pour une personne du metier 

•&■ document qui fait partie de la meme famille de brevets 



Date a laqueOe la recherche Internationale a ete effectlvement achevee 



6 novembre 2003 



Date d'expedition du present rapport de recherche Internationale 

13/11/2003 



Nom et adresse postale de r administration charges de la recherche Internationale 
Office Europeen des Brevets, P.B. 5818 Patentlaan 2 
NL-2280HVROswiIk 
TeL (+31-70) 340-2040. Tx. 31 651 epo nl, 
Fax: (+31-70) 340-3016 



Fonctlonnaire autorlse 



Modesto, C 



Formutaire PCT/1SA/210 (deuxi6ma fertile) QuIUet 1992) 



RAPPORT DE REC 



HE INTERNATIONALE 



Demande^^fc|tionale No 

PCT/fIW/01755 



C.(sulte) DOCUMENTS CONSIDERES COMME PERTINENTS 



Categorie 1 Identification des documents cites, avecje cas echeant, Hndlcationdes passages pertinents 



no. des revendlcationsvisees 



US 5 930 730 A (KIRLIN R LYNN ET AL) 
27 juillet 1999 (1999-07-27) 
colonne 5 -colonne 8 
figures 1-10 

US 2002/022930 Al (DALLEY RICHARD MAPES 
ET AL) 21 fevrier 2002 (2002-02-21) 
cite dans la demande 
page 1, alinea 16 -page 4, alinea 59 
figures 1-6 

US 6 151 555 A (PEPPER RANDOLPH E F ET 
AL) 21 novembre 2000 (2000-11-21) 
colonne 3, Hgne 45 -colonne 5, ligne 38 
figures 1-39 



1-10 



1-10 



1-10 



Formulatre PCT/1SA/210 (suite da la deuxfeme feuQIa) Quillet 1892) 



RAPPORT DE RECHERC 

Rensejgnements relatffs aux mem! 



INTERNATIONALE 

"amines de brevets 



Demand) 

PCT/ 



tonale No 

;/01755 



1 ■ 

Document brevet cit6 


1 — 


Date de 




Membre/s\ de la 1 


Date de 


au rapport de recherche 


1 


publication 




famille de brevet(s) I 


publication 


US 5615171 


A 


25-03-1997 


US 


5432751 A 


11-07-1995 






US 


5251184 A 


05-10-1993 








US 


5153858 A 


06-10-1992 








AT 


129577 T 


15-11-1995 








AT 


190136 T 


15-03-2000 








CA 


2088501 Al 


10-01-1993 








CN 


1070267 A ,B 


24-03-1993 








DE 


69205673 Dl 


30-11-1995 








DE 


69205673 T2 


20-06-1996 








DE 


69230731 Dl 


06-04-2000 








DE 


69230731 T2 


03-08-2000 








DK 


548341 T3 


04-12-1995 








DK 


652447 T3 


24-07-2000 








EP 


0548341 Al 


30-06-1993 








EP 


0652447 A2 


10-05-1995 








ES 


2144480 T3 


16-06-2000 








GR 


3033559 T3 


29-09-2000 








NO 


930843 A 


07-05-1993 








RU 


2107931 CI 


27-03-1998 








WO 


9301508 Al 


21-01-1993 


UJ DloOU/U 


A 
M 


24-10-2000 


US 


6014343 A 


11-01-2000 






CA 


2219216 Al 


30-04-1998 








FR 


2755244 Al 


30-04-1998 








FR 


2760275 Al 


04-09-1998 








GB 


2320968 A ,B 


08-07-1998 








NO 


975004 A 


04-05-1998 


US 4672546 


A 


09-06-1987 


CA 


1249876 Al 


07-02-1989 






EP 


0181216 A2 


14-05-1986 








NO 


854447 A 


09-05-1986 


US 5930730 


A 


27-07-1999 


US 


5563949 A 


08-10-1996 






AU 


709621 B2 


02-09-1999 








AU 


7382996 A 


28-04-1997 








CA 


2204168 Al 


10-04-1997 








CN 


1166207 A 


26-11-1997 








EP 


0796442 Al 


24-09-1997 








NO 


971801 A 


05-06-1997 








MO 


9713166 Al 


10-04-1997 








US 


RE38229 El 


19-08-2003 








AU 


696742 B2 


17-09-1998 








AU 


4133396 A 


03-07-1996 








CA 


2179901 Al 


20-06-1996 








CN 


1138902 A ,B 


25-12-1996 ; 








EG 


20609 A 


30-09-1999 








EP 


0736185 Al 


09-10-1996 








NO 


962731 A 


11-10-1996 








RU 


2144683 CI 


20-01-2000 








WO 


9618915 Al 


20-06-1996 








US 


5838564 A 


17-11-1998 








RU 


2187130 C2 


10-08-2002 


US 2002022930 


Al 


21-02-2002 AU 


5616501 A 


03-09-2001 






CA 


2398854 Al 


30-08-2001 








WO 


0163323 Al 


30-08-2001 








EP 


1257850 Al 


20-11-2002 








JP 


2003524192 T 


12-08-2003 



Formutalre PCT/lSA/210 (annexe f am tiles de brevets) Quillet 1892) 



RAPPORT DE RECHERCHE INTERNATIONALE 

Renselgnements reialifs aux memb^^^amllles do brevets 



Demand <^^^ca 

PCT/H^; 



itlonale No 

1/01755 



Document brevet cite 
au rapport de recherche 



Date de 
publication 



Membre(s) de la 
famille de brevet(s) 



Datede 
publication 



US 2002022930 Al 



NO 



20024030 A 



US 6151555 



21-11-2000 



AU 
CA 
EP 
GB 
NO 
W0 



3866300 A 
2362848 Al 
1198776 Al 
2370878 
20014353 



A 
A 



0054207 Al 



23-08-2002 



28-09-2000 
14-09-2000 
24-04-2002 
10-07-2002 
07-09-2001 
14-09-2000 



FeimulBliB PCT/1SA/210 (annexe families de brevets) <MDet 1892) 



